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Abstract 

We show that QCD contains stable four-quark QQqq hadronic states in the limit where 
the heavy quark mass goes to infinity. (Here Q denotes a heavy quark, q a light antiquark 
and the stability refers only to the strong interactions.) The long range binding potential 
is due to one pion exchange between ground state Qq mesons, and is computed using chiral 
perturbation theory. For the Q = b, this long range potential may be sufficiently attractive 
to produce a weakly bound two meson state. 
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1. Introduction 

QCD with three light quark flavors contains mesons (with qq flavor quantum numbers) 
and baryons (with qqq flavor quantum numbers) in its spectrum. In addition there are 
"nuclei" which are weakly bound states of the baryons. At the present time there is no 
evidence for other types of states that are stable with respect to the strong interactions. 
It was originally suggested by Jaffe that there should be hadronic resonances with 
qqqq flavor quantum numbers and there is evidence that some of the observed hadronic 
resonances should be interpreted in this way. In this paper we study the possibility of stable 
exotic QQqq hadrons, where Q is a heavy quark (i.e., tuq 3> Aqcd) and the stability refers 
only to the strong interactions. It is easy to see that such states exist in the htlq — > oo 
limit. For very heavy quarks, the quark pair QQ can form a small color antitriplet object 
of size (a s {ttlq) tuq) 1 with a binding energy of order a 2 s {mo) tuq. The two heavy quarks 
act as an almost point-like heavy color antitriplet source with a mass of about 2rriQ for 
the two light antiquarks in the QQqq hadron. This results in bound states that have the 
light degrees of freedom in configurations similar to those in the observed and A c states, 
with the QQ pair playing the role of the heavy antiquark. In the heavy quark limit, the 
binding energy of the QQ pair tends to infinity, whereas the energy of the light degrees of 
freedom is of order Aqcd ? so that the QQqq state has a lower energy than two separate 
Qq mesons and is stable with respect to the strong interactions. 

This argument for the existence of exotic QQqq states in the spectrum of QCD is 
based on the short range color Coulomb attraction in the channel 3 <S> 3 — > 3. For the 
case where Q is the top quark, this description is likely to be quantitatively correct and 
establishes the existence of exotic ttqq states that are stable with respect to the strong 
interactions. The charm and bottom quarks are not heavy enough for their short range 
color Coulomb attraction to play an important role in the formation of a QQqq state, since 
a s( m Q) m Q i s n °t large compared with Aqcd- If QQQQ states exist for Q = c or 6, they 
may be weakly bound two meson systems and the formation of such a bound state depends 
on the potential between the lowest lying Qq mesons. At long distances this potential is 
determined by one pion exchange and is calculable in chiral perturbation theory. For the 
remainder of this paper we examine the picture of QQqq hadrons as two weakly bound Qq 
mesons, and apply it to the c and b quark systems. 

Section 2 contains a derivation of the long range potential using chiral perturbation 
theory. In Section 3 the eigenstates of the potential operator are classified and Aqcd/^q 
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corrections to the Hamiltonian are discussed. In Section 4 we present a variational cal- 
culation that suggests there might be a weakly bound state involving B and B* mesons. 
Some concluding remarks are also given. 



2. The Long Range Potential Between Heavy Mesons 

In the limit rriQ — > oo, the angular momentum of the light degrees of freedom, s#, 
is a good quantum number. Mesons containing a single heavy quark come in degenerate 
doublets with total spins s± = si ± 1/2. The ground state multiplets with Qq a flavor 
quantum numbers (qi = u, qi = <i, q% = s) have si = 1/2 and negative parity, giving a 
doublet of pseudoscalar and vector mesons that we denote by Pa and Pa respectively. 
For Q = c these are the (-D , D + , D s ) and (D , D* + , D^) mesons and for Q = b they are 
the (B~ : B : B S ) and (B*~, B*) mesons. 

The interactions of these heavy mesons with the 7r, K and 77 is determined by the 
chiral and heavy quark symmetries of the strong interactions. The effective Lagrangian 
that describes the low momentum strong interactions of the pseudo-Goldstone bosons with 
the Pa Q ^ and Pa^ mesons is 

+ l JL TvH (Q) H^^ 5 [^d^ - £8"?] 

+ Ai Tr H iQ) H {Q) \^m q i + tfm q g] + X[ TrH (Q) H (Q) Tr[m q E + m g E f ] 
+ ^_ TilI (Q) (T ^ H (Q) 

m Q 

where the ellipsis denotes terms with more derivatives, more factors of the light quark 
mass matrix 

(m u \ 
m q = m d , (2.2) 
\ m 8 j 

or more factors of 1/rriQ associated with violation of heavy quark spin symmetry. The 
traces in eq. (|2.1| ) are over flavor and spinor indices. The pseudoscalar and vector heavy 
meson fields Pa, , Pajl are combined to form the 4x4 Lorentz bispinor matrix 

H (Q) = (l+A [P <Q) 1 ,_ P (Q) l5h (2 .3) 



The field 77 (Q) destroys P {Q) and P*^ mesons with four- velocity v M . The subscript v on 
H'Q\ p(Q> and has been omitted to simplify the notation. The conjugate "barred" 

field is defined by 



H 



(Q)c 



The pseudo-Goldstone bosons appear in the Lagrangian through 



.iM 
£ = exp — 



(2.4) 



(2.5) 



where 



M 



K~ 



7T 



K 



and 



(2.6) 



(2.7) 



In eq. (^.5|), / is the pion decay constant, / ~ 132 MeV. 

Under SU{3)l x S77(3).r chiral symmetry transformations 

where L G S77(3)z, and i? G SU(3)r. The transformation law of £ is then 



(2.1 



(2.9) 



where L7" is a complicated function of L, and the mesons M. In general C/ depends 
on space-time, but for transformations in the unbroken SU(3)v subgroup V = L = i?, 
and (7 is the constant matrix V. Under SU(2) V heavy quark spin symmetry 
and 577(3) z, x SU(3)r chiral symmetry, the heavy meson fields transform as 



(2.10) 



where S G SU(2) V is the heavy quark spin transformation. 

The long range potential between heavy mesons is determined by the one pion ex- 
change Feynman diagram in fig. [1], with the virtual pion having momentum transfer 
q v = (Q,g). At low momentum, the coupling of the ?W).p*W) and P*«)-P*(Q) to the 



Goldstone bosons is obtained by expanding the Lagrangian eq. (|2.1| ) up to first order in the 
pion fields. The only term which contributes to the one-pion coupling at low momentum 
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transfer is the term proportional to g in eq. ( |2.1|) . In the ttiq — > oo limit the and 
p*{Q) are degenerate in mass and can be treated as a single li H particle." This gives the 
iif-pion interaction 

Ant = -j Tr H iQ) H^ lul5 d v ix . (2.11) 

— * 

The interaction can be rexpressed in terms of the spin of the light degrees of freedom Si 
and the isopsin I A of the H field as 



C 



hit 



2ygg 
/ 



(^■<9tt a ) I a . (2.12) 



The Fourier transform of the interaction potential obtained from eq. ( |2.12| ) between two 
H particles is 

f _ _ (Sn • qj (§12 ■ q 

Y h ■ h =T5— ; 

j z q z +m 



V, (q) = — 72 h- h v '„2 ■ ( 2 - 13 ) 



7i 2 denotes the isospin of heavy meson 1 and 2 and Sn^ denotes the spin of the light 
degrees of freedom in heavy meson 1 and 2. In coordinate space, eq. (|2.13|) gives the 
potential 



V n (x) = AI X ■ I 2 
where 



Si X -xSi 2 -x- \ S a ■ Si2 ) W 2 (r) + S £1 ■ S e2 W Q (r) , (2.14) 



^w = ^ e " m "(^ + ^ + v) ' (2 - 15) 



and 

The coupling g determines the D* — > Dit decay rate. At tree level 

r(^+^ J D°7r+) = ^|^| 3 . (2.17) 

The decay width for D* + — > D + n° is a factor of two smaller (this follows from isospin 
invariance). The experimental upper limit on the D* + width of 131 keV []3| combined 
with the measured branching ratios for D* + — > D + n° and D* + — > D°n + leads to the 
limit g 2 < 0.5. A measurement of the branching ratio for D* + — > D +/ y could also lead 
to valuable information on g ||]. The axial current obtained from the Lagrange density 
eq. (O is 

q T A 7,75 q = -g TtHHj v j 5 T a + ... , (2.18) 
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where the ellipsis denotes terms containing one or more Goldstone boson fields and T A is 
a flavor SU (3) generator. Treating the quark fields in eq. ( |2.18| ) as constituent quarks and 
using the nonrelativistic constituent quark model to estimate the D* matrix element of the 
l.h.s. of eq. (|2.18| ) gives g = 1. (A similar estimate of the pion nucleon coupling constant 
gives gA = 5/3.) In the chiral quark model there is a constituent quark pion coupling 
f|. Using the measured pion-nucleon coupling to determine the constituent quark pion 
coupling gives g 2 ~ 0.6. Thus, various constituent quark model calculations lead to the 
expectation that g is near unity. 

Equation ( |2.14j ) is the leading contribution to the long distance part of the heavy 
meson interaction potential. Loops and higher derivative operators give contributions 
that are suppressed by factors of (Airfr)~ l . It seems reasonable that eqs. ( j2.14| )- (|2.16| ) 
dominate the potential at distances greater than r min = [l/2m^) ~ 0.7 fm. In the nuclear 
potential the corrections to one pion exchange are important even at this distance 0. 
We believe this is (at least partly) due to integrating out the A resonance which has a 
large coupling to Nn. In the heavy meson case both the an d P*(Q) mesons are kept 
in the lagrangian (|2.1] ). The lightest heavy mesons that are integrated out of the theory 
do not couple very strongly to P^n and P*(Q'n since in the constituent quark model 
they correspond to P-wave orbital excitations. At r = (2m n )~ 1 , W2 (r m ; n ) = 518g 2 MeV 
and Wq (r min ) = 9g 2 MeV, where we have used the neutral pion mass in computing the 
numerical values. 

The eigenvalues of the potential operator eq. ( |2.14[) are easily determined. The po- 
sition space part of the state vector can be taken to be \z), corresponding to the spatial 
wavefunction 8 3 (x — rz), since the potential is rotationally invariant. Then eigenstates of 
V n (x) are \I I 3 )\z)\K k) where 



is the total spin of the light degrees of freedom (the spin quantization axis is also taken to 
be the z-axis) and 



K = Six + S12 , 



(2.19) 




is the total isospin. Acting on these eigenstates the potential energy eq. ( |2.14J ) can be 



rewritten as 
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TABLE 1 



I 


K 


k 













9 Wo/4 





1 










1 


±1 




1 








-3W /4 


1 


1 







1 


1 


±1 


x x 



S\ x S z n is 1/4 for states with k = ±1, and is —1/4 for states with k = 0, so we obtain the 
eigenvalues of the potential for the states as given in Table 1, 
where X and X\ are defined by 



and 



Xt 



W2 

3 



Wo 
4 



m Wo 

6 4 



(2.22) 



(2.23) 



Since Wq, Xq and Xl are positive, the attractive channels have (I, K, \k\) = (0, 1, 1) , (1, 0, 0) 
and (1,1,0). At r min the potential energies for these states are about — 266g 2 MeV, 
-7g 2 MeV and -170# 2 MeV respectively. 



3. Classification of Eigenstates 

In the previous section we found spatial ® spin parts of the eigenstates of the potential 
V that were of the form \z)\Kk). By rotational invariance R (g) [ l-R' fc) ] is an eigenstate 
of V with the same energy for any rotation g. It is convenient to combine these states into 
ones with a definite "angular momentum of the light degrees of freedom" F usingB J7[ 

\Ff;Kk) = V2F + 1 [ dg D% } *(g) R(g) [ \z)\Kk) ] , (3.1) 

JSU{2) 

1 F is the total angular momentum minus the spin of the heavy quarks. It is not the true angu- 
lar momentum of the light degrees of freedom because it contains the orbital angular momentum 
of the heavy quarks. 
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( F) 

where D fk J (g) is the rotation matrix in representation F and the measure is chosen so 

dg = 1. (3.2) 



fk 

that 



>SU(2) 

Alternatively we can combine states with definite orbital angular momentum 

\£m) = VWTl J dg D<j${g) R(g)\z) , 
with the spin of the light degrees of freedom to get states 

\Ff;£S) = J2^r;Ss\Ff) \£r)\Ss). (3.4) 

r,s 

Using eqs. ([0) - (ELD ^ ^ s straightforward to show that 



(3.3) 



9/4-1 

(/■" f *\ /■'/: /v k) = H,-r \r f *k s \/ (F f\£0; K k) 

= 5 F/F S r f5 K s(-l) K+k (£0\K - k-Fk) . 

This allows a partial wave decomposition of the eigenstates of the potential. Consider 
for example states with K = 1. Then S = 1 and so we can form states with orbital angular 
momentum £ = F — 1 , F, F + 1 . The other restriction is that F > \k\, so that the D matrix 
in eq. (|3.1|) exists. For definiteness, consider the case F = 1. Then for k = we have 



according to eq. (|3.5|) the partial wave decomposition 



i*=o> = Vi li=0) ~yl li = 2) ■ (3 - 6) 



For k — ±1 it is convenient to form the linear combinations 



which decompose as 



\k = 1) ± \k = -1) , 
l±> = J ^ 1 , ("J 



l+> = Vl l^=°> + Vl l^ = 2 > ' (3, 



-) = \l=l). 
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TABLE 2 



F 


K = S 


k ± 


l 


K- for / = 1 


V n for I = 














-3 Wo/ 4 


9W /4 


1 








1 


-3W /4 


9W /4 


2 








2 


-3W /4 


9W /4 





1 





1 


-X 


3X 


1 


1 





0,2 


-x 


3X 


1 


1 


1+ 


0,2 




-3Xx 


1 


1 


1- 


1 


X 1 


-3X 1 


2 


1 





1,3 


-X 


3X 


2 


1 


1+ 


1,3 


X! 


-3Xx 


2 


1 


1- 


2 


X! 


-3Xi 



Table 2 gives the eigenstates of the potential and their eigenvalues up to F = 2. 

The "angular momentum of the light degrees of freedom" F must be combined with 
the spin of the heavy quark pair Sq = 0, 1 to get the total angular momentum J of the 
bound state. In addition, if the two heavy quarks are of the same flavor, then only states 
which are completely symmetric are allowed. This gives the additional restriction that 
I + £ + Sq is even. Since Wo, X and X\ are all positive it is easy to identify the channels 
which have an attractive potential. In the tuq — > oo limit, the kinetic energy of the heavy 
mesons can be omitted. In this case, QQqq bound states exist if and only if there exists 
some channel in which the potential is attractive. It is obvious from Table 2 that there 
exist attractive channels, so we have another demonstration that there exist exotic states 
in the heavy quark limit of QCD. 

There are two Aqcd/wq corrections to the Hamiltonian that are important for the 
case of heavy but finite quark masses, such as for the b or c quark. The kinetic energies of 
the heavy mesons 



^kin = 



2^ 



where 



1 



1 



1 



+ 



(3.9) 



(3.10) 



is the reduced mass, should be included. In addition, at order Aqcd/^q the — p*(Q) 
mass difference A^) = m P *(Q) — m P (Q) should be taken into account. Experimentally 
A^ ~ 141 MeV and A( b ) ~ 46 MeV so for Q = c and b this effect is quite significant. 
We define the mass splitting so that it is zero on p( < ?) states and on P*(<2) states. 

This adds a term to the Hamiltonian, A, that is diagonal on the p(Q)p(Q\ p(Q)p*(Q) j 
p*{Q)p*{Q) basis. The relationship between this "meson-type" basis and the \Kk)\SQ sq) 
basis is straightforward to determine. The meson basis is obtained by first combining the 
spins of the heavy and light degrees of freedom in each H particle, whereas the K, Sq 
basis is obtained by first combining the two heavy spins and the two light spins. A 
straightforward computation of this change of basis gives: 

|1,1)|0,0) = I/|p*(<2i),i)|p(Q2)} - |p(Qi))|p*(Q 2 ) 7 i) 

+ \P< Q i\l)\P*( Q *\0) - |P*( Ql \0)|P*^ 2) ,l)} , 

|1,0)|0,0) =-{|P* (c?l) ,0)|P (Q2) ) - |p(Qi))|p*(Q2) 5 ty 

21 (3.11) 
_l_ |p*(Qi) ; i)|p*(<22) ; _^ _ |p*(Qi) ) _i)|p*( < 32) ; ]_^| 

|1,-1)|0,0) =i{|p*(<3i),-i)|p(Q2)) _ \p(Qi))\p*(Qi),-i) 

- |P*( Ql) ,-l)|P* (Q2) ,0) + |P* ((3l) ,0)|P* (c?2) ,-l>} . 

The mass splitting term is not diagonal in the K, Sq basis used in Table 2. The general 
form of the mass splitting in this basis is complicated, and is discussed in the Appendix. 
However, it is simple to compute the expectation value of the mass splitting in a given K 
state. For example, 

(K=l,k = ±1\A\K = 1, k = ±1) = ^ (A( Ql ) + A( Q2 ) + 2A^ Q ^ + 2A^ 2 )) . (3.12) 

We need to obtain the energy "meson molecules" including both the 

kinetic energy and the mass splitting, to determine whether they are bound. A particularly 
promising entry in Table 2 is the / = 0, Sq = state on the sixth line. This has the most 
attractive potential allowed in the table, and has a small orbital angular momentum barrier. 
It is also a state which is allowed when the two heavy quarks are identical. In the next 
section we examine the energy of this state in the case when both heavy quarks are 6-quarks 
using a variational calculation. This state is not an eigenstate of the Hamiltonian, since the 
mass splitting and the kinetic energy are not diagonal in the K basis. The mass splitting 
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will produce mixing to other \Kk) \Sq sq) states, which will only lower the energy further, 
and make the state more strongly bound. This F = 1 state has total angular momentum 
one and even parity so it cannot decay (strongly) to BB. As long as the expectation value 
of the Hamiltonian 7i = Tikin + V + A in this state is less than A^ (the mass of a widely 
separated B and B*) it is stable with respect to the strong interactions. The state can 
decay electromagnetically to BB'j if the expectation value of the Hamiltonian is positive. 
Otherwise, it can only decay by the weak interactions. 

There are Aqcd/i^q effects we have neglected. For example at this order the heavy 
quark symmetry relation between the P*(Q> p*(Q)^ and p(Q> p*(Q)<^ couplings is altered. 
This effect however is less important than those we have included. 
4. Variational Calculation of the Binding Energy for Q = b 

In this section we consider the case where both heavy quarks are 6-quarks and focus 
on the state in line 6 of Table 2 with 1 = and Sq = 0. The energy E (above 2ms) of 
this trial state is 

poo 

E[4>]= / r 2 dr 0* (r)H(j)(r) , (3.13) 



where 

1 



n 







d 2 | 2 d 
dr 2 r dr 



m B 

4> (r) is a trial radial wavefunction normalized to unity. 



+ -^-a +V(r) + l (m* B - m B ) • (3.14) 



oo 

2 



r l dr 0* (r) if) = 1 . (3.15) 



The orbital angular momentum barrier follows from eq. ( |3.ij| ). The state being considered 
is a linear combination of I = and I = 2 with probabilities 2/3 and 1/3, so the effective 
value of L 2 is (0)(0 + l)(2/3) + (2) (2 + l)(l/3) = 2. Thus the angular momentum barrier 
is the same as for an I = 1 state. 

At large distances the potential V (r) is given by the one pion exchange potential 

V v (r) = -3X 1 (r) . (3.16) 

We expect V (r) to be given by eq. ( |3.16| ) for r > r m ; n = (l/2m n ). However, some 
information on the short range part of the potential is needed to make further progress. 
In the case of nuclear forces there is a short range repulsive core that is often attributed 
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to vector meson exchange 0. The situation in the heavy meson case is quite different. 
Suppose the nonet of qq vector mesons V£ h is coupfed to the heavy mesons via the term 



C = g v Tr H a H b v^ b . (3.17) 
This is the type of coupling the constituent quark model suggests is appropriate. The 



coupling in eq. ( 3.17 ) gives rise to the contribution 



^® = ^T<{l( f2 - 3 / 2 ) + i} ■ < 3 ' 18 ' 

to the potential, where we have taken m p = = my, and defined I to be the total 
isospin in the two particle channel. The first term in the braces comes from p exchange 
and the second comes from 10 exchange. Note that u exchange is repulsive as in the case 
of the nuclear potential. However, while in the nucleon case the w-exchange piece is much 
larger than the rho exchange piece, for the heavy meson potential rho exchange dominates 
in the 1 = channel giving an attractive potential from vector meson exchange. There is 
no repulsive hard core in the heavy meson potential. Multiple pion exchange contributions 
to the heavy meson potential are also less important than for nucleons. The pion-nucleon 
coupling constant is g 2 = 1.56, whereas the heavy-meson nucleon coupling is g 2 < 0.5. 
This, and the fact that the analog of the A resonance is not integrated out, implies that 
multiple pion graphs are relevant in the heavy meson case only for much smaller values of 
r than in the nucleon case. 

Given that vector meson exchange as modeled by eq. (|3.18|) is attractive in the channel 



we are considering, a conservative approach is to use in our variational calculation the 
potential 

V{r) = { (3.19) 
14 (r) r > r min , 



corresponding to flattening out the one pion exchange potential in eq. ( |3.16|) for r < r m i n . It 



is important to remember, however, that eq. (|3.19| ) is a (conservative) guess and conclusions 
drawn from it should not be taken too seriously. There are physical effects that increase the 
potential energy which we have neglected. For example, r\ exchange gives the contribution 

2 2 (&i " 9j (Se2 ■ q) 

Vv (0) = -TT72 V ^ 2 J S 2 > ( 3 - 2 °) 
3/ 2 q 2 + m* 
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to the Fourier transform of the potential. In the 1 = channel it has the opposite sign 
from pion exchange. Its effects, however, are quite small since in the 1 = channel since 
it is suppressed numerically by a factor of 1/9. In position space, there is an additional 
suppression factor because the potential falls off exponentially in a distance vti~ x rather 
than m~ l . There are also contributions from derivative vector meson couplings to the 



heavy mesons. These are less important than eq. (|3.17|) at large distances but may be of 
comparable importance at r ~ \jra p . In the case of the nuclear potential their contribution 
to the tensor force is thought to be very important even at distances as large as 1 fm. 

The variational calculation using the energy function eq. ( |3.13| ) with Hamiltonian 
eq. (|3.14|) and potential eq. (|3.19|) is straightforward. We have chosen to do the computation 
with g 2 equal to the present experimental bound of g 2 = 0.5. A simple trial wavefunction 
can be chosen of the form 

0( r ) = N e~ ar r b (1 + cr) , (3.21) 

where iV is a normalization constant chosen so that (p satisfies eq. ( |3.15|) . The minimum 
of the energy is at a = 6.23 m^, b = 2.26 and c = —0.16 m^. For these values of the 
parameters the wavefunction </>(r) is peaked near r = r min . The state is bound, with a 
binding energy of 8.3 MeV relative to the BB* energy. (Recall that the state we are 
considering cannot decay into BB.) The average radial kinetic energy is 25.5 MeV, the 
average angular momentum barrier energy is 31.3 MeV, and the average potential energy 
is —88.2 MeV. The mass splitting contributes an additional 23 MeV of energy, since 
the state is 50% BB* and 50% B* B* . 

The binding energy is sensitive to the precise value of g 2 . For example, with g 2 = 0.6, 
it is 26.9 MeV, whereas for g 2 = 0.4, the state is not bound by about 8.2 MeV. The value of 
the binding energy is also sensitive to the value of the potential below r m i n . We have chosen 
to use a flat potential below r m ; n as a conservative extrapolation, and used the neutral pion 
mass in the numerical computations (which gives a weaker potential). The state could be 
much more strongly bound if the potential is more negative than our estimate. We have 
also investiagated the possibility of DB and DD bound states. With the approximations 
we have made, these states are not bound. The reduced mass of these states is small 
enough that the kinetic energy overwhelms the attraction due to the potential. However, 
it is possible that these states are bound if the interaction potential is more attractive than 
our estimate. 

It is interesting that there is a limit of QCD in which one can show that there must 
exist states with exotic quantum numbers. These QQqq states are very difficult to produce 
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or detect experimentally. It is much easier to produce a meson-antimeson bound state. 
The one-pion exchange potential for meson-antimeson bound states is the negative of the 
potential for meson-meson bound states that we have studied in this paper. The meson- 
antimeson spectrum can be investigated by similar methods to those used here. There is 
one important difference — the meson-antimeson sector has annihilation channels which do 
not exist in the meson-meson sector, so there will be no stable bound states. However, 
there might exist resonances. 

We thank F. Feruglio, N. Isgur, S. Koonin and H.D. Politzer for useful discussions. 
This work was supported in part by the U.S. Dept. of Energy under Contract no. DEAC- 
03-81ER40050 and Grant No. DEFG03-90ER40546 and by a NSF Presidential Young 
Investigator Award PHY-8958081. 

Appendix A. Transformation of Basis 

The transformation between the \Kk) \Sqsq) basis (which gives eigenstates of the 
potential) and the usual angular momentum-meson type basis is computed explicitly in 
this appendix. The states \Kk) are obtained by combining the spins of the light degrees 
of freedom in the two meson, and the states \Q q) are obtained by combining the spins of 
the heavy quarks in the two mesons. (We will denote the states \Sq sq) by \Q q) from now 
on, to avoid multiple subscripts in the formulas.) It is useful to define the states 

\Pp-KQ) = Y,(Kk-Qq\Pp) \Kk)\Qq) , (A.l) 

k,q 

where the angular momentum P is the sum of K and Q. It is convenient to treat the spin 
zero meson p(*2) and the three possible polarizations of the spin one meson P*(Q) using 
a unified notation. For this reason, let |Wi«;i) denote the first meson, where |0 0) is the 
spin zero P^\ and |1 w) denotes the three possible polarization states of the vector P*(*2) 
meson. The two meson state is then denoted by \W\ w\) \W 2 W2). Define the state \S s) to 
be the state obtained by combining the total spins of the two mesons, 

\Ss;W 1 W 2 } = (Wiw 1 ;W 2 w 2 \Ss) \W l w 1 }\W 2 w 2 } . (A.2) 

Wl ,W2 

The transformation formula? will involve the overlap of the two states \Pp;KQ) and 
\S s; Wi W 2 ). Now \Pp; K Q) is obtained by first combining the light spins of the mesons 
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into K and the heavy spins of the mesons into Q, and the resultant into P, whereas 
| S s; W\ W 2 ) is obtained by first combining the light and heavy spins of the first meson 
into W±, and of the second meson into W 2 , and the resultant into S. The overlap is 
therefore 



(Pp;KQ\Ss;W 1 W 2 ) = ^{2W X + 1) (2W 2 + 1) (2K + 1) (2Q + 1) 

( 1/2 1/2 W t ) (A.3) 
x I 1/2 1/2 W 2 } 5 PS 5ps, 
{ K Q S 

using the definition of the 9-j symbol. 

The eigenstates of the potential with total angular momentum j are 

\j m;KkQq) = y^j + l [ dg (g) R(g) [\z)\K k) \Q q) ] . (A.4) 

JSU(2) 

These states are not identical to the ones defined in eq. ( pO| ) because we have also in- 
cluded the heavy quark spin \Q,q) along with \K,k) in the definition of the state. The 
conventional states are obtained by taking the spatial wavefunctions of definite orbital 
angular momentum eq. (|3.3| ) and combining them with the spin state of the mesons given 
by eq. 

\jm;£SW 1 W 2 ) = v / 2/Tl ^(£r;Ss|jm) 

r r (A - 5) 

x / dgD®*(g) R(g)\z) \S s;W 1 W 2 ) . 

JSU{2) 1 J 

The transformation matrix is then 

(jm;£SW 1 W 2 \jm;KkQq) = ^{21 + 1) (2j + 1) ^(lr-Ss\jm) 

, ^ (A.6) 

x / dgD^ k (g)D^(g) (S s;W 1 W 2 \ R(g) \K k) \Q q) . 

JSU(2) 



Substituting the inverse of eq ( |A.1| ) into eq. (|A.6| ), inserting a complete set of states 



\P' p'; K' Q') to the left of the rotation operator, and using 

(P'p';K'Q'\R(g) \P p; K Q) = D$ (g) 5 PP , 8 KK > 5qq> , (A.7) 
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leads to the expression 

(jm;£SW 1 W 2 \jm;KkQq) = 

y/(2£ + l)(2j + l) (Zr;Ss\jm)(Kk;Qq\Pp) 



r,s,p,p' ,P 



(A.8) 



x / dg D ( £l +q (g) D%)D$(g) (S s;W 1 W 2 \Pp';KQ) 

JSU{2) 



Rewriting the products of D matrices using the Clebsch-Gordan decomposition, using the 
orthogonality of the D matrices, and using eq. ( |A.3|) reduces the above expression to 



7^777(1/2 1/2 

(jm;£SW 1 W 2 \jm;KkQq) = J f— ' { 1/2 1/2 



(2j + 1) \ ' K 



Q 




(A.9) 



x y/(2W! + 1) (2W 2 + 1) (2K + 1) (2Q + 1) 
x (Kk;Qq\Sk + q) (10; S k + q\j k + q) . 

The 9-j symbols are invariant under reflection about either diagonal. In addition, the 
symbols are invariant under even permutations of the rows or columns, and are multiplied 
by (— l) s under odd permutations of rows or columns, where E is the sum of all nine 
parameters. Thus the independent 9-j symbols that we need for our problem are 

1 r 1/2 1/21] 1 



1/2 


1/2 


1/2 


1/2 








1/2 


1/2 


1/2 


1/2 








1/2 


1/2 


1/2 


1/2 





1 



2 n i V54 



1 



1/2 


1/2 


1/2 


1/2 





1 


1/2 


1/2 


1/2 


1/2 


1 


1 


1/2 


1/2 


1/2 


1/2 


1 


1 



2>/3' n 1§! 



1 

6' 



(A.10) 



1/2 1/2 1 "J 
1/2 1/2 1 
112 



9' 



Using these values we can compute the decomposition of the various \j m; K k Q q) states, 
e.g. 

|lm;1100) =^ |lm;0110) - ^- |lm;0101) - A- |lm;0111) 

-^|lm;1110) + ^|lm;1101) + i|lm;llll) (A.ll) 
+ ^=ilm;2110)-^=ilm;2101)-^|lm;2111) 
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where the states on the right hand side of the equation are \jm;£SWi W%) states. The 
decomposition of the \j m) state with K = 1, k = — 1, Q = q = can be obtained from 
eq. ([A by changing the sign of all the I = 1 terms. Thus the state K = 1, k = + which 



\lm;l + 0) = |lm;0110) - 4g |lm;0101) - |lm;0111) 

+ -^=|lm;2110)-^=|lm;2101)-^|lm;2111). ^ 

Rewriting the Wi and VF2 labels using the more familiar P^ and P* (Q) labels gives 



From this decomposition, it is easy to see that the state on the left hand side is 25% 
p(Qi)p*(Q 2 ) j 25% p*WOp(Qa) and 50% P*(Q^P*(Q^ : and is 67% £ = and 33% £ = 2. 
The left hand side of eq. ( |A.13|) was expressed in terms of angular momentum eigenstates 
in eq. (|3.8|) , and in terms of meson states in eq. ( |3.11|) , which are special cases of the 
simultaneous decomposition in eq. ( |A.13| ). 



is l/y/2 times the sum of the k = 1 and k = — 1 states can be decomposed as 




(A.13) 



lm;21P^p*«a)^ _ lm;21P*( Ql )p*W a )V 
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Figure Captions 

Fig. 1. The one-pion exchange contribution to the meson- meson potential. The solid 
lines are either the or P*^\ and the dashed line is the pion. 
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